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We study the generation of large scale primordial magnetic fields by a cosmological phase transi¬ 
tion during the radiation dominated era. The setting is a theory of N charged scalar fields coupled 
to an abelian gauge field, that undergoes a phase transition at a critical temperature much larger 
than the electroweak scale. The dynamics after the transition features two distinct stages: a spin- 
odal regime dominated by linear long-wavelength instabilities, and a scaling stage in which the 
non-linearities and backreaction of the scalar fields are dominant. This second stage describes the 
growth of horizon sized domains. We implement a recently introduced formulation to obtain the 
spectrum of magnetic fields that includes the dissipative effects of the plasma. We find that large 
scale magnetogenesis is very efficient during the scaling regime. The ratio between the energy den¬ 
sity on scales larger than L and that in the background radiation r(L,T) = ps{L,T) /p cm t(T) is 
r(L,T) ~ 10 -34 at the Electroweak scale and r(L,T) ~ 10~ 14 at the QCD scale for L ~ 1 Mpc. 

The resulting spectrum is insensitive to the magnetic diffusion length. We conjecture that a similar 
mechanism could be operative after the QCD chiral phase transition. 
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I. INTRODUCTION 

A variety of astrophysical observations including Zeeman splitting, synchrotron emission, Faraday rotation mea¬ 
surements (RM) combined with pulsar dispersion measurements (DM) and polarization measurements suggest the 
presence of large scale magnetic fields [J 1, 0, 3, 0, 3- The strength of typical galactic magnetic fields is measured 
to be ~ /i Gj2| 0, 3 @ an d they are correlated on very large scales up to galactic or even larger reaching to scales 
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of cluster of galaxies ~ 1 MpciSil. The origin of these large scale magnetic fields is still a subject of much 
discussion and controversy. It is currently agreed that a variety of dynamo mechanisms are efficient in amplifying 
seed magnetic fields with typical growth rates T ~ Gyr -1 over time scales ~ 10 — 12 Gyr (for a thorough discussion 
of the mechanisms and models seeQ). The ratio of the energy density of the seed magnetic fields on scales larger 
than L (today) to that in the cosmic background radiation, r(L) = pb(L)/ p cm b must be r{L ~ IMpc) > 10 -34 for a 
dynamo mechanism to amplify it to the observed value, or r(L ~ IMpc) > 10 -8 for the seed to be amplified solely 
by the gravitational collapse of a protogalaxyQ. 

There are also different proposals to explain the origin of the initial seed. Astrophysical batteries rely on gradients 
of the charge density concentration and pressure and their efficiency in producing seeds of the necessary amplitude 
is still very much discussed^, Q|. Primordial magnetic fields that could be the seeds for dynamo amplification can 
be generated at different stages in the history of the early Universe, in particular during inflation, preheating and 
or phase transitions[4, Q. Primordial (hyper) magnetic fields may have important consequences in electroweak 
baryogenesisjjjj, Big Bang nucleosynthesis (seejjl Y. the polarization of the CMB^j via the same physical processes 
as Faraday rotation, and structure format ion [J |5, jjjj, thus sparking an intense program to study the origin and 
consequences of the generation of magnetic fields in the early Universe {jjJ-fM). 

A reliable estimate of the amplitude and correlations of seed magnetic fields must include the dissipative properties 
of the plasma, in particular the conductivity pH flil ITil . In ref. [23] we have introduced a formulation that allows to 
compute the generation of magnetic fields from processes strongly out of equilibrium. This formulation, which is based 
on the exact set of Schwinger-Dyson equations for the transverse photon propagator is manifestly gauge invariant and 
is general for any matter fields and any cosmological background (conformally related to Minkowski space-time). In 
the case in which strongly out of equilibrium effects arise from long-wavelength fluctuations, such as during phase 
transitions, this formulation allows to separate the contribution of the hard degrees of freedom which are in local 
thermodynamic equilibrium from that of the soft degrees of freedom that fall out of LTE (local thermal equilibrium) 
during the phase transition and whose dynamics is strongly out of equilibrium. This separation of degrees of freedom 
leads to a consistent incorporation of the dissipative effects via the conductivity (for details seej^). I n that reference 
a study of magnetogenesis in Minkowski space-time during a supercooled phase transitions was presented and the 
results highlighted the main aspects of the generation of magnetic and electric fields in these situations. 

The goals of this article: In this article we study the generation of large scale (hyper) magnetic fields by 
a cosmological phase transition during a radiation dominated era by implementing the formulation introduced in 
ref. [25| . The setting is a theory of N charged scalar fields coupled to an abelian gauge field (hypercharge). We 
consider the situation when this theory undergoes a phase transition after the reheating stage and before either the 
Electroweak or the QCD phase transition, since we expect that these transitions will lead to new physical phenomena. 
The non-perturbative dynamics out of equilibrium is studied in the limit of a large number N of (hyper) charged 
fields and to leading order in the gauge coupling. The non-equilibrium dynamics of the charged scalar sector features 
two distinct stages. The first one describes the early and intermediate time regime and is dominated by the spinodal 
instabilities which are the hallmark of the process of phase separation and domain formation and growth. This stage 
describes the dynamics between the time at which the phase transition takes place and that at which non-linearities 
become important via the backreaction. The second stage corresponds to a scaling regime which describes the slower 
non-equilibrium evolution of Goldstone bosons and the process of phase orderingj26| and growth of horizon-sized 
domains. This scaling regime is akin to the solution found in the classical evolution of scalar field models with broken 
continuous symmetries after the phase transition that form the basis for models of structure formation based on 
topological defects [27, [28|. 

The solution of the scalar field dynamics (2(| is the input in the expression for the spectrum of the magnetic field 
obtained in [25| to obtain the amplitude of the primordial seed generated during both stages. 

We find that scaling stage is the most important for the generation of large scale magnetic fields. Large scale 
magnetic fields are generated via loop effects from the dynamics of modes that are at the scale of the horizon or 
smaller. The resulting spectrum is rather insensitive to the diffusion length scale which is much smaller than the 
horizon during the radiation dominated era. The ratio of the magnetic energy density on scales larger than L (today) 
to the energy density in the background radiation r(L,?y) = ps(L,ri)/p cm b{v) is summarized in a compact formula 
|ea. H4.M4|l ]. For L ~ 1 Mpc (today) we find r(L,rj) ~ 10 -34 at the Electroweak scale and r(L,rj) ~ 10~ 14 at the QCD 
scale, suggesting the possibility that these primordial seeds could be amplified by dynamo mechanisms to the values 
of the magnetic fields consistent with the observed ones on these scales. 

In section II we introduce the model, in section III we summarize the dynamics in the different stages after the 
phase transition and discuss the dynamics of gauge fields including the dissipative effects of the plasma. In section IV 
we compute the spectrum of the primordial magnetic field generated during the different stages and assess the regime 
of validity of the approximations invoked. Our results and conclusions are summarized in section V. 
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II. MAGNETIC FIELDS IN FRIEDMANN-ROBERTSON-WALKER COSMOLOGY 

The cosmological setting in which we are primarily interested corresponds to a symmetry breaking phase transition 
in a radiation dominated Universe. Such phase transition is in principle different from the electroweak one 1 and 
presumably occurs at a much higher energy scale, such as the GUT scale ~ 10 15 Gev but is assumed to be described 
by a particle physics model that includes many fields with (hyper)-charge either fermionic or bosonic. We will not 
attempt to study a particular gauge theory phenomenologically motivated by some GUT scenario, but will focus our 
study on a generic scalar field model in which the scalar fields carry an abelian (hyper)charge. The simplest realization 
of such model is scalar electrodynamics with N charged scalar fields (f y, r = 1,... ,N and one neutral scalar field ip 
whose expectation value is the order parameter associated with the phase transition. 

This model is inspired by the 0(4) linear sigma model which is the low energy effective theory of QCD that describes 
chiral symmetry breaking and the dynamics of pionsf29j. 130]. In this low energy effective theory the neutral field 
associated with ip ~< qq > acquires an expectation value while the three pion fields 7r , n° are the (quasi)-Goldstone 
modes associated with chiral symmetry breaking. The charged pions couple minimally to the electromagnetic field |2flj 
and obviously chiral symmetry breaking (a nonzero expectation for the neutral field) preserves the gauge symmetry. 
We argue later that the model [see eq. ED] can describe magnetogenesis during the QCD phase transition. The 
mechanisms of magnetogenesis discussed in the present article is therefore akin to the photoproduction during the 
nonequilibrium chiral phase transition [3~il| . 

The neutral field is not coupled to the gauge field and its acquiring an expectation value does not break the 
17(1) gauged symmetry. This guarantees that the abelian gauge symmetry identified with either hypercharge or 
electromagnetism is not spontaneously broken to describe the correct low energy sector with unbroken U(1)em- We 
will take the neutral and the N complex (charged) fields to form a scalar multiplet under an 0(2N + 1) isospin 
symmetry. As the neutral field acquires an expectation value this isospin symmetry is spontaneously broken to 
0(2N). The explicit breaking of the 0(2N + 1) symmetry induced by the electromagnetic coupling further reduce 
this symmetry to SU(N) x 17(1). If the neutral field acquires a non vanishing expectation value, the isospin symmetry 
breaking does not affect the masslessness of the photon (it will obtain a Debye screening mass from medium effects). 

The action that describes this theory in a general cosmological background is given by 

<r d v ip + v,<p* V v ^j + M 2 (Y + (Y + ***) - \ T ^ ^ 

( 2 . 1 ) 

where 

ieA^ and = d^A u - d^A^ . (2.2) 



and 


N N 

tf<P = y ^4>l<pr , D^D^cp = y^(<9 M + ieA^cpl (d^ - ieA^)cp r . 

r =1 r =1 

Furthermore, anticipating a non-perturbative treatment of the non-equilibrium dynamics of the scalar sector in a 
large N expansion, we have rescaled the quartic coupling in such a way as to display the contributions in terms of 
powers of 1 /N, keeping A fixed in the large N limit. 

A Friedmann-Robertson-Walker line element 


ds 2 = dt 2 — a 2 {t) dx 2 , (2.3) 

is conformally related to a Minkowski line element by introducing the conformal time 77 and scale factor C{ji) as 

11 = j a(t) ’ = ’ ( 2 ' 4 ) 

In terms of these the line element and metric are given by 

ds 2 = C 2 {rf) {drj 2 — dx 2 ) ; = C 2 (rj) , (2.5) 


1 


If the electroweak phase transition is weakly first order, nucleation will be almost indistinguishable from spinodal decomposition and 
the phenomena studied here may be of relevance. 
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where r]^ v = diag(l, —1, —1, —1) is the Minkowski metric. Introducing the conformal fields 

x) = x) ; x) = C{r]) x) , 1 < r < N ; ^(?j, x) = C{rj) ip(t(r]),x) 


and in terms of the conformal time, the action now reads 


S = drj d 3 x 


with 


LIU 


-dS&d v * 


■D^*D V $ 


— M 2 {rj) ( 


_ 


47V 


'k 2 \ 1 

— + $*$) - -F I1V F a0 rr r, 


LLU 


M\i]) = -p 2 C 2 ( v ) - 


C"(y) 

C{v) 


Dfj. — ^ii ie.A^ 


F, 


— A A _ f) A 


( 2 . 6 ) 


(2.7) 


and the primes refer to derivatives with respect to conformal time. Obviously the conformal rescaling of the metric 
and fields turned the action into that of a charged scalar field interacting with a gauge field in flat Minkowski space- 
time, but the scalar held acquires a time dependent mass term 2 . In particular, in the absence of electromagnetic 
coupling, the equations of motion for the gauge held A M are those of a free held in hat space time. This is the 
statement that gauge helds are conformally coupled to gravity and no generation of electromagnetic helds can occur 
from gravitational expansion alone without coupling to other helds or breaking the conformal invariance of the gauge 
sector. The generation of electromagnetic helds must arise from a coupling to other helds that are not conformally 
coupled to gravity, or by adding extra terms in the Lagrangian that would break the conformal invariance of the gauge 

helds [HI- 

The conformal electromagnetic helds £, B are related to the physical E , B helds by the conformal rescaling 


E = 


£ 


C 2 (v) ’ 


B = 


B 


C 2 {rj) ’ 


( 2 . 8 ) 


corresponding to helds of scaling dimension two. A gauge invariant formulation leads to the following Lagrangian 
density (for details seej^S 3^) 

C = ^4- d ^ + ^dfjAr ■ d^A T + ^(VA 0 ) 2 - M 2 (r/) ^'k 2 + Q'k 2 + 

-ieA T ■ - V$ f $) - e 2 (A| - A 2 ) - ie A 0 , (2.9) 

where <L is a gauge invariant local held which is non-locally related to the original helds, and At is the transverse 
component of the vector held (V • At = 0) and Aq is a non-propa gati ng held as behts a Lagrange multiplier, its 
dynamics is completely determined by that of the charge density |25l l.‘S2l . 

The main point of this discussion is that the framework to obtain the power spectrum of the generated magnetic 
held presented below is fully gauge invariant. 

The theory described by the Lagrangian eq. 12.611 above bears some similarity to the scalar-gauge held theory that 
describes semilocal strings 33, 34| l35l| , however there are important differences between the model studied here and 
that studied in refs. ItI l34l l35l |: i) we assume that symmetry breaking occurs along the neutral direction thus the 
charged scalar held does not acquire an expectation value, whereas in the semilocal theory of refs. B H H the 
charged helds acquire an expectation value, and the gauge symmetry is spontaneously broken, ii) We study the 
dynamics beginning from an initial state in LTE above the critical temperature, follow the dynamics through the 
phase transition and compute systematically to lowest order in a em the non-equilibrium spectrum of magnetic helds 
generated by the process of phase separation. The goal of the studies in ref. [3 are very different focusing on the rate 
of production of semilocal strings. The initial state studied in these references places the scalar held at the minimum 
of the (classical) potential and the phases are distributed at random, with particular initial conditions on the gauge 
helds, namely conditions corresponding to zero temperature, broken symmetry states. Furthermore the dynamics in 
ref. |34l is studied in hat space time with an ad hoc dissipative term for the scalar held. 

In contrast, we study the full quantum dynamics beginning from a state of LTE above T c evolving the quantum 
Heisenberg equations of motion and calculate the magnetic held consistently to lowest order in a. Thus while the 
theory studied here and that proposed in refs. 13311341 .1311 bear a resemblance, they describe very different physics and 
we study a different set of phenomena. 


2 Here we neglect the effect of the conformal anomaly 1211 
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III. PHASE TRANSITIONS IN RADIATION DOMINATED COSMOLOGY 

A. Kinematics 

We consider a phase transition corresponding to the dynamics of small field models where the scalar field has 
vanishing expectation value but with a symmetry breaking potential, namely at the top of the potential hill. 

In a radiation dominated cosmology, the initial state is that of local thermodynamic equilibrium at an initial 
temperature T >> T c . Using finite temperature field theory in an expanding background geometry, it is shown 
that the effective time dependent mass term depends on the effective time dependent temperature T(t) = T/a{t) 
which reflects the cooling from the cosmological expansion (see below). Hence at a given time the temperature equals 
the critical and the phase transition occurs. Field modes with wavectors much larger than the symmetry breaking 
scale g will remain in LTE and will not be affected by the symmetry breaking dynamics 

We normalize the scaling factor C(r/) at the reheating time rj = r/R in such a way that C(tir) = 1 then the explicit 
expression for C{if) reads 


C(V) =H R r] 


(3.1) 


where Hr is the Hubble constant at the reheating time, Hr = . 

We can relate Hr to the reheating temperature and the Planck mass G~ 5 through the equation 


and the Einstein-Friedman equation 


* 9* rri4: 


1/2 T 2 

H «=[r G ?) =it 


(3.2) 


(3.3) 


where g* is the effective number of degrees of freedom at the reheating temperature and we introduced the scale M* 
of the order of the Planck mass 


M* = 


3 y/5 1 


2 7r 2 y/g* G 


In radiation dominated epoch the time-dependence of the mass term cm I is given by the expression 

—M 2 (r]) = g 2 H\ tj 2 = ry 2 , 
where we see the emergence of a new mass scale 


(3.4) 


(3.5) 




= \J M Hr . 


(3.6) 


This scale will play an important role in the following discussion and in the comparison with results obtained in 
Minkowski space-time[25|. There is a last scale which plays a relevant role, the horizon scale ruirf) which is fixed by 
the evolution on the time of the Hubble constant: 

rH i.v) = = C{v) V = H R V 2 ■ (3-7) 

Modes with physical wavelength \ p h ys = k 2 ™ inside the horizon 

Xphysig) ~ k~ kys {g) <t h {v) (3.8) 

are causally connected, modes outside the horizon are causally disconnected. 

The relaxation rate of hard modes of the charged fields is given bvpul] 

r(?y) ~ a T(g) In — , (3.9) 

a 

where the effective temperature varies with time as 


T(v) 


Tr 

C(g) 


(3.10) 
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and the expansion rate given by 


Therefore, 


H( Jl) 

rfa) 

H(v) 


T 2 (y) 

M* ' 

10 16 

T(? ? )[Gev] ' 


(3.11) 

(3.12) 


Thus hard modes are in thermal equilibrium for T{rf) < 10 15 Gev. 

In particular, modes with k ~ Tr are the hard modes that give the leading contribution to the conductivity in 
the high temperture limit [3^, (3^- Modes with k <a will manifest the long-wavelength spinodal instabilities and 
their dynamics will be strongly out of equilibrium iJ E 0 0 . Their amplitude becomes non-perturbatively 
large ly, till EH E3 and will be responsible for the non-equilibrium generation of the primordial magnetic field [25$. 

Using eqs. ED and emu we can write the conformal time as 


Tr 

V Hr T(v) ' 

As it will become clear below an important cosmological quantity is the product 


The ratio 


k r] 


k 

CM 


C(v) V 


kphysjv) 

T(V) 


T(v) r H {v) = 


2ir M* 

ATr Tfa) ' 


kphysiyi) _ 27T 

T(v) ~ 


(3.13) 


(3.14) 


(3.15) 


is a kinematical invariant. Its value today is determined by the scale L which will be typically chosen to correspond 
to a galactic scale or the scale of galaxy clusters, and the temperature of the CMB. It is given by, 


LTk = 3.7 x loK (db) . (3.HS) 

Therefore, 

10 -22 for T(rj) = T R ~ 10 15 Gev 

10 -9 for the EW transition (3-17) 

10 -6 for the QCD transition 


; T( V ) V L ) 


for L ~ 1 Mpc. Thus, during the regime of interest in this article, kq 1 for scales of galaxy clusters. A noteworthy 
aspect of eq. (irrm is that the wavelengths corresponding to the scale of galaxies or clusters today were well outside 
the horizon during the radiation dominated era when the electroweak and QCD phase transitions occurred. 

Another important quantity is the ratio of the wavevector k of the primordial magnetic field to the conductivity. 
As it will be discussed below, the physical conductivity is given by 


cr(ry) 


C N( V ) T( V ) 
a ln ^y 


(3.18) 


where C is a constant of 0{ 1), N(r)) is the number of ultrarelativistic charged species, and we have neglected the 
(logarithmic) dependence on the energy scale in the running coupling constant. For this discussion we will neglect 
the time dependence of N(rf) assuming that the number of charged ultrarelativistic species remains constant (this 
assumption can be relaxed without qualitative modifications of the main argument). Under this assumption 


r(? ? ) = 


C(V) 


with an being the comoving conductivity determined at the time of reheating 

CNTr 


vr = 


aln Sv 


(3.19) 


(3.20) 



















Thus the ratio, 
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kphysjll) 27TQ 27 ( M P C \ 

<xfa) ~ NLT r ~ \ L J 

neglecting logarithmic corrections. 

Furthermore, 


°(v) r H(v) 


M. 

aT{r]) 


10 5 for T{rj) = T R ~ 10 15 Gev 
10 18 for the EW phase transition 
10 21 for the QCD phase transition 


(3.21) 


(3.22) 


where we have neglected logarithmic corrections. Therefore ctr ij 1 throughout the radiation dominated era 
considered in this article. The regime <jr 77 3> 1 ; kr iff ctr -C 1 is dominated by the (slow) hydrodynamic relaxation 
of the magnetic field. 

Another relevant estimate involves the (comoving) diffusion length £,diff(v ) = \Zv/ a R 


Zdiffjy) ~ l uT(ri) 

v y m* 


10 " 3 for T{rf) =T r ~ 10 15 Gev 
10 -9 for the EW phase transition 
10 _1 ° for the QCD phase transition 


(3.23) 


where again we have neglected logarithmic terms. Therefore the diffusion length is much smaller than the Hubble 
radius during the radiation dominated era. Finally, combining eqs. iron and (ETT7I) we find 

10- 25 <fce dl// (77)<10- 16 , (3.24) 


between reheating and the time of the QCD phase transition. 

The contribution from the hard modes of both the charged scalar and gauge fields which remain in local ther¬ 
modynamic equilibrium lead to an effective mass for the scalar field. This thermal mass is obtained from the long- 
wavelength limit of the scalar field self-energy and includes the hard thermal loop contributions from the gauge and 
scalar fields S3, B] ■ This thermal mass is given by 

™ 2 t = || (A + 3 e 2 ) . (3.25) 

Finally, another important quantity is the Debye screening length that determines the scale at which long-range forces 
are screened by the polarizability of the medium. In an ultrarelativistic plasma, the comoving Debye screening length 
is given bv fi^B 

(3.26) 

e T r 

the ratio of the Debye screening length to the Hubble radius is given by 


1 T(rj) 
du e M* 


(3.27) 


Hence £d <C d//(? 7 ) for T[gf) < 10 16 Gev, thus long range forces are screened over very short distances. The formation 
of long-wavelength domains with typical size of the order of the Hubble radius [26J leading to strong charge and current 
fluctuations that will seed magnetic fields, will not be hindered by long-range forces, which are effectively screened 
over sub-horizon distances. 

Magnetic field generation via charge asymmetries during a period in which electromagnetism was spontaneously 
broken was previously studied by Dolgov and SilkBl w ho argued that long-range forces would be screened by the 
Higgs mechanism. This is different from the situation studied in this article, where the U (1) symmetry associated with 
electromagnetism (rather hvpercharge) is not spontaneously broken. Long range forces are screened by the plasma, a 
situation not considered in 45|. 


B. Scalar fields dynamics 

For completeness and to highlight the aspects of the non-equilibrium dynamics most relevant to the generation of 
magnetic fields, we summarize the main features of scalar field dynamics. For further details the reader is referred 












to [51Ell!?il- In what follows we will neglect the backreaction of the gauge fields on the dynamics of the scalar fields. 
The rationale for this is that the main non-equilibrium processes that lead to magnetogenesis will be non-perturbative 
in the scalar sector and result from the instabilities associated with the phase transition. The contribution from 
the gauge fields, in the form of self-energies for the scalar fields, do not feature the instabilities associated with the 
phase transition and will, furthermore, be suppressed at least by one power of a the (hyper) electromagnetic coupling 
constant as compared to the scalar self-interaction. 

As described above, the non-equilibrium evolution of long-wavelength modes begins with the spinodal instabilities 
which result in an exponential growth of the amplitudes for long-wavelength fluctuations. When the non-linearity 
becomes of the same order as the tree-level terms in the equations of motion, the back reaction of these fluctuations 
shuts-off the instabilities US, El 112 • Therefore a non-perturbative treatment of the dynamics is required. The large 
N limit of the scalar sector allows a systematic non-perturbative treatment of the dynamics which is renormalizable 
and maintains the conservation laws|40ll4dll. 

We will therefore study the dynamics in leading order in the large N limit that already reveals the important 
non-equilibrium features of the evolution. 

Radiative corrections: 

The contribution from the gauge fields to the equations of motion of the long-wavelength modes of the scalar fields 
arise through self-energy corrections. To lowest order in a these are dominated by the hard modes of the gau ge fi elds 
with momenta ~ T (hard thermal loops) which lead to a contribution to the thermal mass given by eT / y/8 |43| |44| . 
Thus the lowest order radiative corrections had already been accounted for in the thermal mass eq. Id. 2511 . 

The non-equilibrium effects in the gauge contribution of the scalar self-energy will arise from polarization loops in 
the photon propagator. This in turn will induce non-equilibrium radiative corrections to the self-energy of the scalar 
fields of the order 0(a 2 ). These small contributions can be safely neglected in this context. Thus, to this order the 
radiative corrections to the scalar field from the gauge field propagator in the scalar self-energy are accounted for in 
the thermal mass. 

Hence the dynamics of the scalar field is studied along the same lines as presented in refs. minni but the only 
difference is in the initial conditions in the modes that reflect the thermal mass in LTE. 

Since symmetry breaking is chosen along the direction of the neutral field T we write 

^(x, rj) = VN ip(rj) + x(x, rj) ; {x(x, 77 )) = 0 (3.28) 

where the expectation value is taken in the time evolved density matrix or initial state. The leading order in the large 
N limit is obtained either by introducing an auxiliary field and establishing the saddle point or equivalently by the 
factorizations |4pl EJ 

(<f > t $) 2 —> 2 ($ t <i>)<i> t $ 




This factorization that yields the leading contribution in the large N limit makes the Lagrangian for the scalar fields 
quadratic (in the absence of the gauge coupling) at the expense of a self-consistent condition: thus charged fields 4> 
acquire a self-consistent time dependent mass. 

The dynamics is determined by the Heisenberg equations of motion of the neutral field and the charged fields 
$ |4(j, l4lL WM . We will consider that at the onset of the radiation dominated era, the system is in the symmetric high 
temperature phase in local thermal equilibrium with a vanishing expectation value for the scalar fields. Consequently, 
the initial conditions are < H/(af, 0) >= y/N y>(0) = 0, < W(x, 0) >= y/N ip{ 0) = 0, < $ r (0,£) >= 0, < 4> r (0,i ? ) >= 

0 . 

In the absence of explicit symmetry breaking perturbations the expectation value of the scalar field will remain 
zero throughout the evolution, thus ip = 0 . 

It is convenient to introduce the mode expansion of the charged fields 


'Mfoz) = 


d 3 k 




Or(k)f k (v ) e ik - s + bt(k) r k ( V ) e- ik -■ 


r = 1,... ,N , 


(3.29) 


with < a r (k) >=< b r (k ) >= 0 . 

In leading order in the large N limit, the Heisenberg equations of motion for the charged fields translate into 
the following equations of motion for the mode functions and the expectation value of the neutral field for 77 > 

riR; mm ei 
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dr/ 2 




<p(v) = 0 • 


(3.30) 


Obviously, the initial conditions <^(0) = 0, (/>(0) = 0 imply that tp(ri ) = 0 for all times. That is, = 0 is a fixed 
point of the dynamics. 

We must now append initial conditions for the mode functions fk(v)- The initial conditions on the mode functions 
fkijl ) depend on the value of the wavevector k as compared to the horizon scale H R l at the reheating time: 

• k > Hr: for fluctuations inside the horizon, we may assume thermal quasi-particle boundary conditions at the 
reheating temperature: 

Mvr) = y= 5 fkiVR.) = ~iW k fk{VR ) , W k = \Jk 2 + m 2 T , (3.31) 

where the frequencies Wk are quasi-particle frequencies with thermal mass m R given by eq. (TT25l at the reheating 
temperature, 


Wk = \Jk 2 + rrij, . 

For these modes the assumption of local thermodynamic equilibrium is well motivated and we have 

(4(fc) a s {k)) = (bl(k) b s {k)) = S rs n k ; n k = —^- (3.32) 

e~ - 1 

• k < Hr: for superhorizon fluctuations, which are causally disconnected at the reheating time, we cannot assume 
a thermalized distribution. The correct distribution has to be derived by following the dynamics from the 
inflationary stage, when the fluctuations were well inside the horizon. While a complete discussion of the initial 
conditions is left to a forthcoming article, the case under consideration we will see that the dependence on the 
initial conditions is rather weak and only during the initial stages of the phase transition. For the later stages, 
dominated by the scaling solution described below, the dynamics is universal and does not depend on the initial 
conditions. We will simply assume that both Wk and nk have a finite non-zero limit as k —> 0 namely the only 
important quantities for the dynamics of long-wavelength fluctuations are 

lim Wk = Wo , 0 < Wq < oo, lim nk = no , 0 < no < °o • (3.33) 

fc—»0 fc—>0 

With this choice of the initial state we find the backreaction term to be given by 

'*> = i /||f lAMI’P + - (3.34) 

This expectation value is ultraviolet divergent, it features quadratic and logarithmic divergences in terms of an upper 
momentum cutoff. The quadratic divergence and part of the logarithmic divergence (the one proportional to the 
mass term) are absorbed in a renormalization of the mass term fi 2 —> fj? R and the remainder logarithmic divergence is 
absorbed into a renormalization of the scalar coupling A —> Ar. While these aspects are not relevant for the discussion 
here, they are mentioned for completeness, the reader is referred to for details. 

After renormalization the self-consistent field is subtracted twice, and is given by (for details see |2(| and 

references therein) 


2 N x ' 

A R I(V) 


A r [J{ri) +I(rj)) , Xr J ( r] ) = ^2 J Q I 2 d( l\fg(v)\ Q 


°q 5 


A R 


S? L q dq ) l/a(,)l “ a + 


1 e(q - K 2 ) 


2 q 3 


~^ 2 r + 2 


(3.35) 


and the mass and coupling are replaced by their renormalized counterparts /j,^; Ar respectively. Here K is an arbitrary 
renormalization scale. In order to avoid cluttering of notation we now drop the subscript R for renormalized quantities, 
in what follows /i; A stand for the renormalized quantities. 

The finite temperature term J(j]) has contributions from short wavelengths for which the mode functions are of 
the form f q (j]) ~ e lqri /\fq and contributions from long wavelengths. The contribution from short wavelengths is the 
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same as that in equilibrium in Minkowski space time and determines the hard-thermal loop l43l El contribution to 
the self-energy given by|H 

Jhtl = (3.36) 

where we have used that the short wavelength modes are in thermal equilibrium at the reheating temperature T R . 
This hard thermal loop contribution has been self-consistently accounted for in the thermal mass of the scalar field 

ea. H3.25D . 

It is convenient to separate the hard thermal loop component eq. iTHol) from eq. lTT35l) and define 


AS fa) 


A 

8-7T 2 



l/?fa )| 2 ( 1 + 2 ««) 


1 

q 



2 

- 1 . 


Q(q-K 2 ) 


2q 3 


-H 2 + A Ufa) 


(3.37) 


After renormalization and in terms of dimensionless quantities, the non-equilibrium dynamics of the charged scalar 
fields is completely determined by the following equations of motion J 2 (| (d(J, [dl], [l 2 j| , 


drf“ 


+ Ad 2 fa) + q 2 + AEfa) 


fa fa) = 0 ; f q (r/ R ) = 




; f'aim) = -iw q fafafa 


with W q given by eq. 13.5D and the effective, (conformal) time dependent mass is given by 


Ad 2 fa) = C 2 fa) /r 2 


Tl 


c 2 fa) fa 2 


- 1 


rj~iZ _ 


24 n 2 
A + 3e 2 ' 


(3.38) 


(3.39) 

(3.40) 


Where we have used <p = 0 as a fixed point of the dynamics. 

The time dependent mass term A4 2 fa) includes the high temperature corrections and clearly displays the cooling 
associated with the expansion in the form of a time dependent effective temperature T e f f fa) = T R /C{rf). The 
phase transition occurs at a time rj c when Tfa^fafa = T c , thus for 77 > rj c the effective time dependent mass term is 
Ad 2 ( 77 ) = M 2 (rj) = —/r 2 H\ rf = —ft , 4 rf as given by equations 13.5D - 13.6D . 

The full time evolution of mode functions in a radiation dominated cosmology has been studied analytically and 
numerically in detail in ref. [26L l4fl| . Here we hig hlig ht the most important features which are necessary ingredients to 
study magnetogenesis. The reader is referred to|2f| for a more comprehensive discussion. 

The are two main dynamical stages in the evolution: 


• Spinodal stage: this is the stage immediately after the phase transition which is dominated by spinodal decom¬ 
position and the growth of long-wavelength fluctuations!^. This stage spans the time scale r] c < rj < rj n i where 
the non-linear time scale rj n i is determined by (see below) 


vh = 


A Efa„z) 


(3.41) 


During this stage the back-reaction, determined by the term A Efa), can be neglected and the dynamics is 
linear. 


• Scaling stage: This is a stage in which the non-linearity encoded by the back-reaction term AS ( 77 ) are very 
important and compete with the tree level term in the equations of motion. This stage is described by a scaling 
solution of the equations of motion for the modes with small wavevectors and describes the non-equilibrium 
relaxation of long-wavelength fluctuations El l27l E3|. 


1. Spinodal stage 


After the phase transition but before the non-linear time scale after which the back-reaction becomes important, 
namely for 77 c <C 77 < r] n i the time dependent mass term is given by Ad 2 = —ft 4 i] 2 , and for weak coupling A< 1 we 
can neglect the back-reaction AEfa). The equations of motion for the mode functions during this stage are given by 



~4 

-M 1 


2 


fq(rf) = 0 ; q < ft 2 v 


(3.42) 
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We note that for T c ~ 10 15 GeV 


A* Vc 


( /A+ 3 e 2 M* 
IV 24 T c 


10 


(3.43) 


and therefore for 77 > ? 7 C we are in the regime ft rj 1 . It is clear that the mode functions f q (rj) will increase 
exponentially in the band of unstable wavevectors q < ft 2 Tj. Eq. iEH5l) can be solved exactly in terms of Hermite 
functions 1401 


fgijl) =b q e " (/i?7)2 H i^_^(^v) + a q e = (/iI?)2 


(3.44) 


where the constants a q and b q are fixed by the initial conditions (|3.38|) . For (i i) > 1 we can use the asymptotic 
behavior of the Hermite functions 491 , 


H v (z) z = l (2 zf 


1 + 0 


1 


and we find for the mode functions, 


f q (v) a q e^ v)2 (ft rj) 


1 1 
— -nr ~ I 


1 + 0 


~ 9 2 

/i Z Vi. 


(3.45) 


Since the exponentially damped solution becomes negligible the phases of the mode functions f q (r)) freeze , namely, 
they become constant in time and are slowly varying functions of q for long wavelengths. 

This is very similar to the situation in Minkowski space-time, where the mode functions however increase as e Mt , 
i.e. much slower. In any case the soft (q —> 0) modes are the most amplified at the end of the evolution, therefore, 
the quantum fluctuations (EH7I) are dominated by the lower integration bound q = 0 . 

We notice that the freezing of the long-wavelength mode functions will play an important role in the discussion 
about the magnetic field generation, since it assures the independence of the final result from the initial particle 
distribution function, except for subleading corrections. 

The physics of the phase transition is essentially the same as in Minkowski space-time|2^, E3 EH E2i si nce the 
exponential growth of modes in the spinodally unstable band will make the back reaction term AE ( 77 ) begin to grow and 
eventually cancel the term —ftftrj 2 in the equations of motion (for 77 >> r] c the effective time dependent temperature 
vanishes). 

This will happen at a non-linear time scale defined by EoL ITHl 

A E (rjni) = ft 4 rj 2 nl (3.46) 


Two important aspects are described by ij n i: i) at this time scale the phase transition is almost complete since 
AE(77 „z) = ft 4 rj 2 a means that \(&+}/2N = ft 4 r] 2 i, namely the mean square root fluctuations in the scalar field probe 
the manifold of minima of the potential. 

ii) At 77 ~ r/ n i the mean square root fluctuations of the field are of order M 2 (r] n i) /A probing the vacuum manifold, 
and the non-linearities become very important. The back reaction AE (?7 n i) becomes comparable to M 2 (i) n i) and the 
instabilities shut-off. Thus for r/ c < r] < rj n i the dynamics is described by the linear spinodal instabilities while for 
V > r] n i a full non-linear treatment of the evolution is required. As it will be discussed below this later stage is 
described by the emergence of a scaling solution. 

For r] n i > rj ft~ 4 the asymptotic form for the mode functions apply and we find for the the quantum 

fluctuations (TQ7I) which dominated by the lower integration bound q = 0 , 


A H(ij n i) — A (1 + 2 no) 


fi 2 |a 0 | 2 


~ 2 2 


1 + 0 


32 7 r 2 ft^ [In (ft rj n i)}~ 

This leads to the following estimate for the spinodal time for weak coupling A 


1 


/i rjnl J J 


Vnl ~2 

+ 2 


In 


32 7T2 


A (1 + 2n 0 ) |a 0 | 2 ft 


In In ( —-————-— ] + O (In In In — 

l A (1 + 2?r 0 ) |a 0 | 2 ji ) \ A 


(3.47) 


(3.48) 
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The important point is that the dependence on boundary conditions and the initial distribution is solely logarithmic , 
thus we may expect out predictions to be very robust with respect to changes of the initial conditions. In particular, 
the scale factor at this non-linear time scale is given by 


= h) 1 

where we have used eqs. and (TH0l) . 

The amplitude of the long-wavelength modes at the non-linear time, roughly speaking at the end of the phase 
transition is approximately 


'"x 


1 + 0 




(3.49) 


\fq(Vnl)\ = 


/ 32 wi 

1 32 7T2 

1 A (1 + 2 77-0 ) fi 

A (1 + 2 710 ) 


1_ <2 
4 2 p. 2 


(3.50) 


As we will discuss in detail below this non-perturbative scale will ultimately determine the strength of the magnetic 
fields generated during the phase transition. 

During the intermediate time regime the equal times correlation function is approximately 


<*£»**,*(»?)> = kl 2 e ^ 2 e-(^ +1 ) lnM . 


and its Fourier transform for long wavelenghts is of the form 


~ 2 2 

(4> a (x, r?)4>b(0, r))) = 6 a , b |a 0 | 2 ky— fi 3 


jlr) U n k 

which determines the time dependent correlation length of the scalar field, 


e z' 2 (v) 


(3.51) 


(3.52) 


£(v) = v / ikk = 2 ' 


I In (y^ H r i f) 
M h r 


(3.53) 


This expression is valid in the intermediate time regime r) c < 77 < r] n i during which the non-equilibrium dynamics is 
dominated by the spinodal instabilities. The detailed analysis of the dynamics in refs. |4(i Dll 142 ] and the discussion 
of the main features presented above can be summarized as follows: 

• At intermediate times /I -1 -C 77 < rj n i ~ / 2 -1 yln 1 /A the mode functions grow exponentially for modes in the 
spinodally unstable band q < M (77 ). The phase of these mode functions freezes , namely, becomes independent 
of time and slowly varying with momentum. 

• At a time scale determined by the spinodal time the back-reaction shuts off the instabilities and the phase 
transition is almost complete. This can be understood from the following: the backreaction becomes comparable 
with the tree-level term (for 77 > r/n) when ik(d>f<I>) « ffif 2 . This relation determines that the mean square 
root fluctuation of the scalar field probes the minima of the tree level potential. 

• During the spinodal stage the correlation length of the scalar field grows in time and is given by eq. 13.5311 . This 
is interpreted as the formation of correlated domains that grow in time, and is the hallmark of the process of 
phase separation and ordering. This correlation length will be important in the analysis of the correlation of 
magnetic fields later. 

• The large fluctuations associated with the growth of spinodally unstable modes of the charged fields will lead to 
current fluctuations which in turn will lead to the generation magnetic fields. Thus the most important aspect 
of the non-equilibrium dynamics of the charged fields during the phase transition is that large fluctuations of 
the charged fields associated with the spinodal instabilities will lead to the generation of magnetic fields. Since 
the modes with longer wavelength are the most unstable the magnetic field generated through the process of 
phase separation will be of long wavelength. Furthermore we expect that the magnetic field generated by these 
non-equilibrium processes will be correlated on length scales of the same order as that of the charged field 
above. 
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2. Scaling stage 


A remarkable result of the evolution in the asymptotic regime (when the effective temperature has vanished) found 
in ref. [26j is that there is a very precise cancellation between the tree level term —]i 2 C 2 {rj) and the back reaction 
AE(r?) in the equations of motion (13.381) . The self-consistency condition requires that for a radiation dominated 
cosmology j26j 


AE^-^C 2 ^) ^ 


15 
4 p 2 


In this asymptotic regime the solutions of the equations of motion 


dr] 2 


+ k 2 


15 

4yy2 


fk(v ) =0 


(3.54) 


(3.55) 


are given by 


fk(v) = Vv 



+ B k k 2 N 2 {kr,) 


(3.56) 


This solution can be written in terms of the scaling variable 

x = k T] 


(3.57) 


in a more illuminating form 


fk(v) = A k v 512 



+ B k 


x 2 N 2 {x) 

jj3/2 ’ 


(3.58) 


As discussed in detail in ref. \2(^ . the relevant integrals are dominated by x ~ 1, namely by modes with wavelength 
of the order of the Hubble radius, thus the second contribution proportional to N 2 (x) can be safely neglected at long 
times. 

For x < 1 in the long time regime we can further approximate A k ~ Ao and the asymptotic solution during this 
stage is of the scaling form 


fk(v) = A q if 12 (3.59) 

x z 

Since for x < 3 and large time the modes with small wavevector have the largest amplitudes, these dominate the 
backreaction. The very precise cancellation it leads to the following sum rule [2l| 


J k 2 dk |/ fc ( 7?)| 2 (1 + 2 n 0 ) ^ =°° ff rf 


(3.60) 


Since for large r] the integral is dominated by soft modes k ~ - —> 0 the distribution function can be approximated 
by no and the amplitude by |A 0 | 2 . The sum rule ea. H3.0Ull then leads to the identity 


where we 


used the integral phi 


IA>| 2 (1 + 2n 0 ) 


30t r 3 


,, 2 FT 2 
]i tl R 



(3.61) 


which is dominated by x < 3. This is a remarkable result: the product \Aq\ 2 [1 + 2 no] in the scaling regime does 
not depend on the initial conditions on the evolution , namely it is universal in the sense that it is independent of the 
previous history through the phase transition. This is an important result which will play an important role in the 
power spectrum of the magnetic fields. 

An important consequence of this scaling solution is that the equal time two-point correlation function of the scalar 
field is given by 


($ a (f,ry)4> 6 (0, p)) = 5 a ,b D(z) 



(3.62) 
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which reveals that the correlation length is given by the size of the causal horizon 26]. The dynamical evolution during 
the scaling stage is precisely determined by the growth of horizon-sized domains [26j. 

We summarize below the important features of the solutions in the scaling regime that will be used in the compu¬ 
tation of the power spectrum of the magnetic fields. 


• For 77 >> r] n i a scaling regime emerges in which the mode functions are given by eq. (Mi with x = k rj. 
This scaling solution describes the relaxation of long-wavelength fluctuations of the charged fields. Again the 
phase of these modes freezes namely is independent of time. This is important because this fact will entail that 
the retarded self-energy of the transverse photon polarization tensor will give a subleading contribution to the 
generation of magnetic fields. 

• The sum rule eq. 13.61)1) constrains the product of the amplitude times the occupation of the long wavelength 
scaling modes to be given by eq. CUB- 

• The scaling solution described above is akin to that found in classical models of formation of topological 
defectsj23,|28|. The scaling regime describes the evolution of long-wavelength fluctuations and the adjustment 
of the spatial correlation length of the scalar field to the Hubble radius |2(|. 


C. Gauge field dynamics 


In a high temperature plasma a very important aspect that must be taken into account in the dynamics of gauge 
fields is the electric conductivity, which leads to dissipative processes. As discussed in j 2 f|, the electric conductivity 
severely hinders magnetogenesis, and also introduces the diffusion length scale which could limit the correlation of 
the magnetic fields that are generated. 

In Minkowski space in equilibrium the conductivity is obtained from the imaginary part of the photon polarization 
and it is dominated by charged particles of momenta p ~ T in the loop with exchange of photons of momenta 
eT < k <C T |3^, [Hj. A careful analysis including Debye (electric! and dynamical (magnetic) screening via Landau 
damping leads to the conclusion that the conductivity is given by jH, EH 


a = 


C NT 

01 ln TJTv 


(3.63) 


with N the number of charged fields and C ~ 0(1). 

In an expanding cosmology an in particular during phase transitions, a more precise assessment of the contributions 
and meaning of the conductivity must be provided. As it was discussed in ref. [ 2 ^], the fluctuations of the charged 
fields during a phase transition will have very different behavior if the typical wavevector of these modes is of the 
order of or smaller than the symmetry breaking scale or much larger than this scale. 

Short wavelength modes, those with typical wave vectors much larger than the symmetry breaking scale are insen¬ 
sitive to the phase transition and are always in local thermodynamic equilibrium (LTE). For short wavelength modes 
deep inside the horizon, the mode functions are of the free field type f q (r]) ~ e lqri /^/q. 

Long wavelength modes, those with wavectors of the order of or smaller than the symmetry breaking scale undergo 
critical slowing down and fall out of equilibrium during the phase transition. These modes become spinodally unstable 
during the early stages of the transition as summarized above and analyzed in detail in refs. 1261 l4fl l41L 143 . 

Thus the contributions from the charged particle fluctuations to the photon polarization must be separated into two 
very different regimes: a) the hard momenta p > \M(rf)\ correspond to charge fluctuations that are always in local 
thermodynamic equilibrium, b) the soft momenta p <C \M(rj)\ fall out of equilibrium and undergo long-wavelength 
spinodal instabilities and enter the scaling regime. 

The contribution from hard momenta will lead to a large equilibrium conductivity in the medium, while the 
contribution to the polarization from soft momenta will contain all the non-equilibrium dynamics that lead to the 
generation of electromagnetic field fluctuations. 

As the instabilities during the phase transition develop, the fluctuations of the charged fields will generate non¬ 
equilibrium fluctuations in the long-wavelength components of the electric and magnetic fields with the ensuing 
generation of long-wavelength magnetic fields. However the large conductivity of the medium will hinder the generation 
of electromagnetic fluctuations, hence the conductivity must be fully taken into account to assess the spectrum of the 
magnetic and electric fields generated during the non-equilibrium staeel25l. 

We are interested in the generation of long wavelength magnetic fields, namely k « T but also k « a 2 T, 
since within the astrophysical application, the wavelength of interest for magnetic fields are of galactic scale, while 
T corresponds to a wavelength at the peak of the CMB ~ millimeters. Thus the physical situation corresponds to 
studying the photon polarization tensor for long-wavelength. 
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The polarization tensor has local (tadpole) and non-local contributions. The equilibrium contribution to the tadpole 
~ (<j)t<j)) is dominated by momenta ~ Tr leading to the hard-thermal loop result ($1*1)) ~ Tj|. This contribution 
is actually cancelled by the zero frequency-momentum limit of the non-local polarization (bubble diagram). This 
cancellation in equilibrium is a consequence of the Ward identities and implies the vanishing of the magnetic mass in 
absence of symmetry breaking |44L 147^ . In equilibrium the inverse propagator vanishes in the limit of zero frequency 
and momentum [441. l47l| in absence of symmetry breaking at any temperature. 

Out of equilibrium, the tadpole and the non-local polarization (bubble diagram) exactly cancel each other in the 
longwavelength, longtime limit as shown in ref. 0 [see eq.(6.12) in the reference] indicating the vanishing of the 
photon mass. 

In equilibrium the long-wavelength and low frequency limit (k, to —> 0) of the spatial and temporal Fourier transform 
of the transverse polarization is given by 

IFr [k, ui) = i u> cr (3.64) 


Thus we write for the full transverse polarization for long-wavelength electromagnetic fields 

IIt(?7, 7/, k) = a- y') + n noneq(ri, 77 ', k) (3.65) 

The first term above includes the contribution from the hard momentum modes p ~ Tr in the transverse polarization, 
while n none q(? 7 , p 1 , k) is the contribution from the long wavelength modes which are unstable in the spinodal stage 
and take the scaling form in the scaling regime. Thus in a very well defined sense, the polarization describes 

the effective low energy theory for the transverse photon field. 

Our strategy is to obtain the non-equilibrium contribution to the spectrum of electromagnetic fields to lowest order 
in a but treating the conductivity exactly. 

In a cosmological space-time, the temperature scales with the inverse of the conformal factor 

nr,) = (3.66) 


and therefore the conductivity a = < 7 ( 77 ) becomes time-dependent. If we are interested in time scales where the number 
of ultrarelativistic charge carriers does not change significantly, which is the case that we will consider in what follows, 
then the time evolution of the conductivity is purely kinematic: 


n„) = m • (3.67) 

An important effect of the conductivity, as discussed in jUJ, is the introduction of a diffusion scale in the transverse 
photon propagator. The long-time behavior of the zeroth order propagators for the transverse gauge fields: retarded 
(R), advanced (A), symmetric (H) 


• v R,A,H( r h’n'r k ) ='P ij (p) T>R,A,H(y,y',k) 

obeyfllj (see[2^ for details) 

^2 + 

^2 + k2 + 

-^ 2 +k 2 + a(i 7) C(p) 

with the transverse projector 

p) = - pV . 


d 

drj 

d_ 

dp 

d 

dr] 


T> R (r), 7?', k) = 5(r] - 7 /) ; T>r(i], 7 /) =0 for 77 < 7 / 
V A (r],T)',k) = 6(r] - 7 ?') ; V A (r], 77 ') = 0 for 77 > 7 / 
V H {r],r)\k) = 0 , 


(3.68) 

(3.69) 


Due to eq. iTToTl) the comoving conductivity or = cr(j]) C(rf) is an invariant quantity in the regime in which the 
number of ultrarelativistic charge carriers is constant. The estimate given by eq. (l3~22l) clearly indicates that during 
the radiation dominated era between reheating and the QCD phase transition, <jr 77 1. 

Then for k <C ctr (which is certainly fulfilled since the relevant wavevectors are k <^T <^T/a ~ a R ) and 77 3> 1/ctr 
we can safely neglect the second order time derivatives in eqs. EHHl. leading to the following equations, 


^R{y,y',k) = V C (77,77'; k) ^(t?-?/) , T> A [r ] 1 r]', k) = T> c {r ] 1 r]'; k) 9 {r]' — 77) 

e -^(nU) e -£(u-i?') 

V H (i],r]',k) =i - , X>c(t 7,?7 ';fc) = 


°R 


&R 


(3.70) 

(3.71) 
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IV. MAGNETIC FIELD SPECTRUM 


As discussed in detail in reference ( 2 f|, the quantity of astrophysical relevance is the correlation function 

<B i ( V ,x)B i (r ] ,0)> p , (4.1) 

where the sum on repeated indices is understood. B(ipx) above is a Heisenberg operator and the expectation value 
is in the initial density matrix. From this quantity, the spectrum of the magnetic field is obtained in the coincidence 
limit 

SB{v,k) = i lim [ d 3 x < {B l (ri,x),B l (r]' , 0)} > p e lk ' x , (4.2) 

Z v ^V J 


where { , } denotes the anti commutator. And from 5 (b(? 7 , &) we can extract the physical magnetic energy density 
stored on comoving length scales larger than a given L 

27T 

Ap B (L,r)) = J k 2 S B {r],k) dk . (4.3) 

where we have restored the powers of the scale factor arising from the transformation to conformal time. Denoting by 
A pb(L,t]) the contribution from the non-equilibrium generation (subtracting the local thermodynamic equilibrium 
contribution), a quantity of cosmological relevance to assess the relative strength of the generated magnetic field is 
given by the ratio of the power on scales larger than L to the energy density in the radiation background 


where 


= Ap B (L,7?) 

Pj(v) 


Pi 
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(4.4) 


(4.5) 


is the comoving energy density in the thermal equilibrium background of photons. 

The physical energy densities ApB, P hys(L , if) , p 1:P h y s are obtained from the comoving expressions above by rescaling 
p —> p/C A {rf) as can be seen from the conformal rescaling 12.811 . Thus the ratio r{L,rj) would be a constant in the 
absence of non-equilibrium generation or dissipative processes. Hence the time dependence of the ratio m only is 
solely a consequence of the non-equilibrium generation mechanisms or dissipative processes (such as magnetic diffusion 
in a conducting plasma) but not through the cosmological expansion. 

Using the results obtained in reference]^! that lead to a first principle derivation of the spectrum, we just quote 
its expression to leading order in a, (here and in what follows Sb refers solely to the non-equilibrium contribution to 
the spectrum) 

Sb(v, k) = Sg(r],k) + Sff {p, k) (4.6) 


with 


S I B (p,k) = e 2 N j cos 2 9) (1 + n q ){l + n |ffjf| ) J dr) X k V c {f), r)u k)f q {vi)f w+ ^{vi) 


+(1 + n q ) 


D n \q+k\ 
V 

VR 


+ n q n \$+k\ 


f dip k r D c {v,m,k)f q (m)f^ + £,(»7i) 

^ VR 


z(l 


n„ l + n 


\q+k\> 


'VR 


dip k T> c (?7,m,fc)/*(77i)/ | - +S| (7 ? i) 


J drp kV c {v,Vi,k)f*(rii)f* s+ k\( r h) 


(4.7) 


and 


S§(v,k) = —ik 2 F(rj, rj; k) , 

where F(r),rj';k) satisfies the homogeneous differential equation 


(4.8) 


j -2 + k 2 + a{rj)C{rj)^- 
drj z dr] 


F(r), r)', k) + J drp [rf ( 771 ) 6 ( 77 - 771 ) +Il R (r],rp)] F(rp,rf,k) = 0 , (4.9) 
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with n^r/), 11 ^( 77 , 77 ') being the one loop tadpole (local) and retarded (non-local) contributions transverse 

polarizationp^|. 

We note that the function F(r), 77 '; k) obeys the same equation as the transverse gauge mean fieldj^, but as it will 
be argued in detail below, its contribution to the spectrum of primordial magnetic fields generated during the phase 
transitions is subleading in the scalar coupling constant A. The equation of motion mu can be solved systematically 
in an expansion in powers of the non-equilibrium polarization, 


£(77,77; k) = F^ (77,77; k ) + F ( ' 1 \ 77,77; k ) + 0(a 2 ) , 

F [0 \ v , v ', k ) = £>#(77,77' ,k) , 

F w {r],r]',k) = f d ‘ n 1 V c ('n ~ Vi, &) [ drj 2 11 *(772) <5(771 - 772) + II# (771,772; k) £>#(772 - 77', fc) + 

JriR Jt]r 

+ (v <->■ v'), 


(4.10) 


where II* , II# are the tadpole (local) and the retarded contribution from the one-loop transverse photon polarization 
respectively (for details seepffh 
These are given by [25| 


II* ( 77 ) = —ie 2 N 


d 3 q 

( 2^)3 


^>(77,77,7) , n#(? 7 i,772, k ) = n > (? 7 i, 7 ? 2 ,fc) - n<(771,772, k ) 0(771-772) ( 4 . 11 ) 


with 


n>( 77 i, 772 ,fc) = 2 ie 2 N 


d 3 q 2 


(2tt) 3 


q (1-cos 6 >) G>( 77 i, 772 ,g) G>(? 7 i, 772 ,|g + fc|) , n<( 77 i,? 72 ,A:) = n>(?72, 77 i,A:) . 


The scalar propagator G>(? 7 , 77 '; k) is constructed from the mode functions f q (r/) that satisfy the mode equations (Id. 4211 
and is given by 


G>(77i,7?2;fc) = - [(1 + rife)/ fc ( 771 )/^( 772 ) +7i fe /fc(?7i)/ fe (?72)] 


(4.12) 


Therefore 


n#(?7i, 772 ,fc) = 4e 2 fV 


d 3 q 

(2^F 


q 2 (1 - cos 2 0) Im G > (r) 1 , t? 2 , q) G > (r] 1 , t? 2 , \q + fc|) 8 ( 771 - 772 ). 


(4.13) 


This expression for the retarded self-energy must be contrasted with that of the contribution from S^(r],k) which 
requires the real part Re G>( 771 , 772 , < 7 ) G > (rji, 772 , \q + k\) . This is an important difference, the long wavelength 

modes of largest amplitude in either phase given by mu' or by (13.5*111 are such that their phases are frozen, namely 
they do not depend on time, therefore the products / 9 (?7i)/|p-+^| ( 772 ) with only the growing mode solutions are real 
and such products will contribute only to S I B ( 77 , k). This freezing of phases is a consequence of the classicalization of 
the scalar field fluctuations 1261. 


We now argue that the contribution from S B is subleading. First of all, the term F^°\r],r];k) in ea. 14.Kill is the 
solution of the homogeneous equation in absence of non-equilibrium fluctuations and leads to the local thermodynamic 
equilibrium contribution to the power spectrum, which is independent of the non-equilibrium generation mechanisms. 
This contribution has been analyzed in detail in ref. j25l| and will be subtracted. In what follows we focus solely on 
the contribution from the non-equilibrium fluctuations. 

For intermediate times after the phase transition during the spinodal stage r] c < 77 < r] n i , the long-wavelength mode 
functions are approximately given by eq. 13.4511 . 

Near the end of the phase transition for 77 ~ 77 ' ~ 77 n i the leading order time dependence of the scalar Green’s 
functions is approximately given by 

fk{ri)fkW) oc i . (4.14) 

where we used eq. mm> - Thus, the contribution from the tadpole (local term in the self-energy) is of the order 

e 2 

n tad ( 77 ) oc — + subleading . 

A 


(4.15) 
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This estimate is consistent with the fact that the tadpole contribution is e 2 < > and near the end of the phase 

transition the mean square root fluctuations of the scalar field probe the vacuum state, namely < |4 >| 2 >~ /i 2 /A. Since 
the phases of these modes are frozen, there is no contribution from the leading order to the retarded polarization, since 
it requires the imaginary part of the product of propagators as displayed in eq. m> . Because of this cancellation 
of the leading term, the contribution from the retarded polarization bubble is of the same order as that of the 


tadpole 12 . 1 . 1 1 


IIr oc y . (4.16) 

A similar argument based on the sum rule C5B leads to the same conclusion in the scaling regime. 

For late times^the k —> 0 limit of the retarded polarization exactly cancels the contribution from the tadpole [see 
eq.(6.12) in ref.|4Sj], 

The contribution from Sj ? is in both cases of 0(1/ A 2 ) since each long-wavelength mode function is of order 1/V\ 
at the end of the spinodal stage or, by the sum rule iTTflOl) in the scaling regime. Thus we can safely neglect the 
contribution from S R to the magnetic spectrum. 


Thus the leading contribution to the power spectrum generated by non-equilibrium fluctuations is given by 


Sb(j 1-> k) = (1 + 2n 0 )' 


aN k 2 2 + r 


e "r 


q A dq d (cos 9) (1 — cos 2 9) 


r^n ,2 


'VR 


e CT « fq(vi) f\ ?+ k\(vi) dm 


(4.17) 


where 9 is the angle between the vectors q and k and where we have replaced 


(1 + 2n q )(l + 2 ri|- + £|) — (1 + 2 n 0 ) 2 , 

since as highlighted in section (ITTTTm the dynamics during both the spinodal stage as well as the scaling stage is 
dominated by the long-wavelength modes that acquire non-perturbatively large amplitudes. 

The final form of the power spectrum generated by the non-equilibrium dynamics given by eq. iITTTI is the basis for 
the study of primordial magnetogenesis during the different stages after the phase transition. 


A. Magnetogenesis during the spinodal stage 


The long-wavelength mode functions in the spinodally unstable band are given by the expression Id.4511 . 

The integral over m f° r large 77 can be computed integrating by parts in eq. ITT71 as an expansion in 1 /(/i. 77 ) 2 . 
The integral is dominated by the upper limit, which leads to the cancellation of the exponentials that contain the 
conductivity. 

The integrals over momenta and angles in eq. (1X171) can be done straightforwardly when the mode functions are 
given by eq. (1X751) . Thus from eq. iIXTTIi whe obtain the following expression for the spectrum of magnetic fields 
generated by the non-equilibrium fluctuations 


S B {k,rj ~ rjni) = 


512 772 N a k 2 
A 2 <j 2 r Z 5 (r] n i) 


0 -\kH 2 (vm) 


1 + 0 


ln I 


where £(??) is given by ea. (15.551) . In obtaining this result we used the following 

[1 + 2 n 0 ] 2 |a 0 | 4 2(mtw) 2 _ 1024 tt 5 
(WnO® ^ A2 d 2 


[see eq. (IXX1 1 and the identities|46l|. 

r +1 


4 q k x 

dx (1 — x ) e + 


in fir] ^ ( 4 q k 

16 (q k In /j. 77) 3 \ fl 2 


ln(/t7y) cosh 


4 q k 


In fli 7 


— sinh 


4gfc _ 

—- in 

d 2 


(4.18) 


(4.19) 


(4.20) 


and 



[£ 2 q k cosh (^ 2 qk) — sinh (£ 2 q k )] 


AA k 3 £e 


\k 2 e 


(4.21) 
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Notice that the magnetic field spectrum tna is independent on the amplitude |ao| and on the initial occupation 
(1 + 2no) 2 - Therefore this result is quite robust. 

This result is the same as for the Minkowski space-time (see eq.(7.47) in ref.@|)i except for a multiplicative factor 
p 6 £ 6 (?? n z) — 8 In 3 (in and the expression for the correlation length in the radiation dominated universe Id. 5411 . 

As in Minkowki space-time, the presence of a high conductivity plasma severely hinders the generation of magnetic 
fields. However, a noteworthy aspect is that up to the non-linear time th e magn etic field is still correlated over the 
size of the scalar field domains rather than the diffusion length £*// ~ \Ji i/<Jr- The diffusion scale determines the 
spatial size of the region in which magnetic fields are correlated in the absence of non-equilibrium generation. The 
ratio between the domain size £(ry) given by eo. lld.5dH and the diffusion length scale £,diff{v) is given by 

t(vm) _ 2 j a R \n(flr] nl ) ~ / M 2 pl In { \ 1 ^ 1 ^ 

tdiffivm) ~ p\ Vnl V M 2 ) 


Where we have used the relations ild.dli . (TOI) . m and Thus an important conclusion of this study is that 

the magnetic fields generated via spinodal decomposition are correlated over regions comparable to the size of scalar 
field domains which are much larger than the diffusion scale. 

The spectrum for the electric field can be obtained from that of the magnetic field by simply replacing k T>c — > 'Dc- 
In the soft regime and for time scales ^ <C we have T>c — —k 2 /cr 2 R whereas kD c ~ k/cr R . Therefore the 
electric field spectrum is suppressed by a factor k 2 /a R with respect to the magnetic field, namely 

S%*(t,k) = ^ S%*(t,k) . (4.23) 

a R 

Thus, in a high temperature plasma with large conductivity the non-equilibrium processes favor the generation of 
magnetic photons instead of electric photons, and again equipartition is not fulfilled. 

The energy density on large scales > L again can be computed in closed form in the limits L £(t n i) and 
L < £(t n i). We find in the first case from ea. 14.3ll . 


Ap B (r]ni,L) 


2 13 tt J # Na 1 

5 A 2 [p ^(r/ni)] 4 o\ £(??„*) L L 4 


T ^ • 


(4.24) 


We find for the opposite case, 


ApsiVnuL) 


3 x 2 10 7 T 3 N a 

AV4FW 


L <C £(t n i) ■ 


(4.25) 


The ratio of the magnetic energy density on scales larger than L at the spinodal time and the total radiation energy 
given by the Stefan-Boltzman law p 1 = 7r 2 T 4 ; /15 results, 


r(r)ni, L) = 


Ap B (VnU L) 3 x 2 13 7T 2 


Na 


P'1 


A 2 [P £(??nz)] 4 <j\ £(Vnl) L (L Tr) 4 


L £(t n l) . 


(4.26) 


This result is the same as for the Minkowski space-time (see eq.(7.53) in ref.[25[), except for a multiplicative factor 
p 6 £, 6 (ri n i) — 8 In 3 (in A) and the expression for the correlation length in the radiation dominated universe 13.5311 . 

The factor (LTr) -4 is purely dimensional and is ultimately the determinining factor for the strength of the generated 
magnetic fields on a given scale. These combinations are invariant under cosmological expansion and are determined 
by the ratio of the scales of interest today (galactic) to the thermal wavelength (today) of the cosmic microwave 
background radiation at the Wien peak. In particular LT R ~ 10 25 for L ~ 1 Mpc(today) [see ea. l,3.1fi[l ]. 

It is clear that the production during this regime is extremely small, due to the large values of (LT) 4 and of the 
ratio Or/ p?. In order to obtain an estimate for the amplitude of the seed magnetic field, we consider the following 
set of parameters: A = 10~ 2 , a = 10 -2 , p = 10 14 GeV, Tr = 10 16 GeV (corresponding to a critical temperature 
T c = 10 1 ' 5 GeV). We then obtain, 


r(L = 1 Mpc) ~ Hr 157 . (4.27) 

Therefore, the amplitude of the magnetic field generated during the spinodal stage is completely negligible. This 
result is similar to the result obtained in Minkowski space-time in ref. |2J| and is expected on the basis of dimensional 
analysis. 
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B. Magnetogenesis from the scaling regime 


In the scaling regime 77 » rj n i the spectrum of the magnetic field is given by eq. iITTtI) with the mode functions in 
the scaling regime given by eq. 111.5911 . 

The final expression for the leading contribution, given by eq. iITTTI reveals a noteworthy aspect. As we have argued 
above, the modes k of astrophysical relevance today, were well outside the horizon during the radiation dominated 
era between reheating and the QCD phase transition. The mode functions eq. (ETT7I) attain the largest amplitude at 
long times for x = qy < 2 — 3, thus momenta in the polarization loop that are within the horizon lead to generation 
of magnetic fields with long-wavelengths well outside the horizon. This we believe, is an important mechanism, loop 
corrections lead to a coupling between modes inside the horizon with those outside. Thus in this manner, causal 
fluctuations can actually lead to the generation of fields with wavelengths much larger than the horizon. 

Since ky -C 1 the power spectrum eq. (1TT7I) takes the following form using the scaling mode functions eq. ESI- 

aN k 2 P P 

S B (y,k) = (l + 2n 0 ) 2 - — \A 0 \ 2 dq / J 2 (qyi) m M / MqM m dy 2 /(<?, k, r?i, y 2 ) , (4.28) 

71 a R Jo J VR JvR 

where we set the exponentials equal to unity in eq. (TT~T7I since ky <C 1 and k <C <J B and 


/(g,fc,77i,?7 2 ) = J dx 


1 — x 2 


-l (<Z 2 + k 2 — 2 kqx) 2 
Using the summation theorem [46| 


J 2 (\/V + k 2 — 2kqx 771 ) J 2 {\/q 2 + k 2 — 2 kqx 772 ) . 


__ A( n 2 1 1, 2 _ Ob nT \ pL 

MVd 2 + k 2 - 2kqxri) = - 2 fc2 ? 2 -^(1 + 2) Ji +2 (qr]) Ji +2 (k 77 ) Cf{x) 

Q ** 1— n 


where the C 2 {x) are Gegenbauer polynomials. For kr] <C 1 the l = 0 terms dominate and we can use the small 
argument behaviour of the Bessel functions J 2 {krj) = ^(A; r/) 2 [1 + 0(k 2 1 7 2 )]. We finally obtain, 

I(q, k, 771 , 772 ) = J 2 (qyi) J 2 O 7772 ) [l + 0(k 2 7? 2 )] . (4.29) 

Inserting eo. 14.2911 into eo. 14.2811 yields 

s B (v,k) = (i + 277 0 ) 2 aN \ \a 0 \ 2 I ^ {? 7 2 [J 2 {.qq) ^ MqV Mqv)] - (?7 -► m)} 2 [1 + e>(fc 2 ? 7 2 )] , (4.30) 

o tt a r j 0 q 

where we used the formula 0] 

z j 2 ((3z) dz = v - [j 2 (Pv) - MPy) MPv)\ ■ (4.3i) 

Since 77 3> 77 ^ we can neglect the terms with rju and we find, for k <C 77 -1 

rv 

S B (r],k) =V -j- — VHrV [1 + 0(fc 2 7? 2 )] ; V = 48.61... . (4.32) 

A a R 

where we used eq. (Turn and we computed numerically the integral 



-5- \Jl{x) - JiOr)J 3 (:r)] 2 = 0.0005295 ... 

X 4 


This integral is dominated by the region x > 1, namely, by modes that are inside the horizon. From the estimate 
E23, the corrections O [k 2 tVnff{rj)) are truly negligible between reheating and the QCD phase transition. 


The dependence on the conformal time ~ 77 ' is a direct consequence of the scaling form of the solution for the mode 
functions. The strong time dependence is a consequence of the causal relaxation of the Goldstone fields, a result of 
the phase ordering kinetics that entails that the size of the domains grow with the horizon. 

The spectrum eq. (1X321) exhibits the following important features: 
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• i: The exponential associated with the diffusion length cancels out, a reflection that the long time behavior of 
the integrals above are dominated by the upper limit. Hence the final result for the spectrum does not feature 
the exponential suppression with the diffusion length. 

• ii: The result for the spectrum only depends on the initial amplitud Aq and initial occupation number no in 
the combination |Ao| 2 (l + 2no) which is constrained by the sum rule en. lid, til 111 . Hence the final spectrum is 
insensitive to the initial conditions on the mode functions or occupations, which in principle carry information 
of the early history beginning from the inflationary stage. This is a consequence of the scaling solution being a 
fixed point of the dynamics of the scalar field 12ft 1271 l28l. 

• iii: A noteworthy result is that superhorizon magnetic fields are generated by the non-equilibrium dynamics of 
modes inside but near the Hubble radius. This is a consequence of the polarization loop, wherein the propagators 
correspond to momenta q and | q + k |. The momenta k corresponding the wavevector (scale) of the magnetic 
field is such that the wavelength is larger than the Hubble radius, but the momenta q corresponding to the 
charged scalar field fluctuations are inside the horizon. The correlation length of the charged scalar field is of 
the order of the Hubble radius. Thus acausal, superhorizon magnetic fields are generated by loop effects. 


In order to reveal the enhancement during the scaling regime in a more transparent manner, it is convenient to use 
the relations (TOli . (THiH) . lEHTili and the explicit expression for the conductivity crnnt in the form 


a R = c[a, N] 


NT r 

a 


c[a,N] — 


ln [alv] 


0 ( 1 ). 


Then, the ratio r(L, ?y) for L rj is given by 


r(r],L) 


2407T T> a 3 
N c 2 [a, N] [LTrY 


VXT( V ) 



(4.33) 


(4.34) 


where T> is given in eq. <l~i~32l) . We note that in the final result there is no dependence on the reheating 

temperature but only on the scale of symmetry breaking /i, the temperature at the time rj and the scalar and gauge 
couplings. This is expected since the non-equilibrium processes begin in earnest after the phase transition, local 
thermal equilibrium prevailed between the time of reheating and the phase transition. 

The dependence on the scalar self coupling oc 1/A 2 is a hallmark of the non-perturbative nature of the growth of 
unstable modes and spinodal decomposition, it is ubiquitous in the non-equilibrium dynamics of phase transitions(2(| 

ini eh 

Large scale magnetogenesis is more efficient for large symmetry breaking scale /r, since the larger the symmetry 
breaking scale, the longer lasts the scaling stage. 

Consider for instance the case in which the symmetry breaking scale /i ~ 10 13 Gev and A ~ a ~ 10~ 2 , corresponding 
to a critical temperature of order of a GUT scale T c ~ 10 15 GeV and suppose that the scaling regime lasts until the 
electroweak phase transition scale, i.e. r/ is such that T(rj) = T EW ~ 10 2 GeV. Then the factor 


/ m. y 

v^A Tew J \Tew ) 


10 


100 


compensates for the factor ( LT R ) 5 . Taking N and <?* of the order of 10 (these values are taken as representative and 
they can be changed simply in the final expressions) we can write the expression for the ratio as 


r(T(i7),L)~10- 34 


1 Mpc 


-5 



Therefore 


r(T(v),L) 


10 34 at the EW transition 
10 -14 at the QCD transition . 


(4.35) 


(4.36) 


Thus if the scaling regime lasts until a time between the EW and the QCD phase transitions the amplitude of the 
large scale magnetic fields is within the range necessary to be amplified by some dynamo models. The amplitude of 
the seed magnetic field is strongly dependent on the duration of the scaling regime. We have only focused on a scaling 
regime terminating either at the EW or QCD phase transition since there will surely be new phenomena associated 
with these that must be included in the dynamics of magnetogenesis. 











C. Discussion 


Validity of the approximations: There are two main approximations that were used to obtain the results 
quoted above, i) the long-wavelength approximation kr) <C 1 and ii) the weak coupling approximation. We now 
provide an estimate of the reliability of both these approximations to establish the limit of validity of our results. 

i): Long-wavelength approximation: In order to reach our final result for the rate r(L, rf) we have explicitly 
used a series of approximations which are valid for long wavelengths but whose validity must be checked before 
we reach any conclusion regarding the spectrum at small scales. In particular we must address the limits of 
applicability of the result eq. (IP1) . This result has been obtained by integrating the magnetic spectrum on 
scales 0 < k < k m ax with k ma x = 27 t / L m i n \ the formula for the magnetic spectrum was valid in the limit 

fomax T]max < 1 . (4.37) 

In order to provide an estimate may take for i) m ax to be the (conformal) time at which the EW phase transition 
occurs, namely t)ew ~ 1 GeV -1 . As discussed in the introduction, we are considering a situation in which the 
magnetic field is considered as a perturbation of a pre-existing thermal blackbody background. For consistency 
this requires that, 


n ^Y J min ; Vwe) 1 . 


(4.38) 


CNa 3 ( T c \ 4 (TewY 
24 2 c 2 \Tew) V ) 


1 -1/5 


This relation translates in a condition 

kmax ‘Hmax 

For T c ~ 10 15 GeV this gives k ma x Vmax 0.0176 which in turns translate into 


A » L r 


70 fm . 


(4.39) 


(4.40) 


However, this is the comoving length normalized at the reheating time. In order to convert to the present time, 
we have to take in account the redshift 


z R = ^ ~ 4 x 10 28 for T r ~ 10 15 Gev ; 
To 


(4.41) 


this gives 


Lmin\today ~ 0.1 P C ' (4.42) 

Thus, the approximations invoked are reliable to estimate the amplitude of primordial seeds on galactic scales 
or larger, today. 

ii): Weak electromagnetic coupling: In order to study the amplitude for much smaller scales the calculations 
must be done without the long wavelength approximations invoked above. In this case we must expect a 
breakdown of perturbation theory and we cannot give a reliable estimate in the present framework. Furthermore, 
for scales well inside the Hubble radius, microphysical processes not included in our approximations, such as 
scattering between charged fields and between charged and gauge fields must be included. These processes 
will restore equilibrium between the different fields, if there is a substantial transfer of power from the charged 
fluctuations to the radiation field, this may lead to a change in the equation of state and the full backreaction 
on the metric must be included. At longer time scales the effects of the backreaction of the gauge fields on the 
dynamics of the scalar field, as well as the non-equilibrium contributions to equation of state and the Friedmann 
equations must be included self-consistently. 

Generation on short distance scales: For scales well inside the horizon during the scaling regime, namely 
qri 1, we must account for causal microscopic processes that tend to equilibrate the electromagnetic fields 
generated by the non-equilibrium processes. In order to understand these processes we must look at the kinetics 
of equilibration. The mode functions for wavevectors well inside the horizon are Minkowski-like, of the form 
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where the coefficients a q , P q must be determined from a full numerical evolution. However the constancy of the 
Wronskian entails that 


Kl 2 - \Pq\ 2 = 1 

which suggests the identification \a q \ = 1 + A f q ; \P q \ = Af qi Af q is the number of (asymptotic) quanta created 
during the time evolution. This form of the asymptotic mode functions leads to the equipartition between the 
electric and magnetic field generation, since spatial and time derivatives are the same. In turn this entails that 
we can understand the generation of electric and magnetic fields by obtaining a kinetic equation for the number 
of photons. Such kinetic equation must necessarily be of the form 

= [! + N k (n)] r>(„) - N k ( V )T<(r,) 

which displays the familiar gain minus loss contributions in terms of the forward and inverse rates. Eventually 
a steady state will be reached which will describe a stationary distribution of photons. The computation of 
the forward [T^(r?)] and inverse [T) 4 ^)] require a detailed knowledge of the distribution A/" g [48| since these 
generalized rates are functionals of these occupation numbers. Clearly such computation lies beyond the scope 
of this article and is a task that we will undertake elsewhere. However, the kinetic equations above will tend to 
an equilibrated state of local thermodynamic equilibrium. 

• Effect on the LSS: It is important to estimate the effect of the magnetic field on scales corresponding to those 
of the last scattering surface, which today are Llss ~ 100 Mpc. From eq. (107)1) we see that at the electroweak 
temperature t(Tew, Llss) ~ 10 -44 , taking the fourth root we can provide an estimate of the temperature 
fluctuation induced by the primordial magnetic field ST/T | igg ~ [t(Tew, Less)]* ~ 10~ n which is negligible 
compared to the CMB temperature fluctuation at this scale ~ 10 -5 . On the other hand, a similar estimate at 
the time of the QCD phase transition gives dT/T\ LSS ~ 10~ 6 which is marginally compatible with the current 
observations. Thus the reliability of the approximation of weak gauge coupling combined with the effects on 
the temperature anisotropy at the last scattering surface seem to lead us to conclude that if a phase transition 
during a radiation dominated era occurs near the GUT scale and results in a scaling stage, our results for 
primordial magnetogenesis will be reliable down to the scale of electroweak symmetry breaking. 


V. CONCLUSIONS 

In this article we studied large scale primordial magnetogenesis during a phase transition in the radiation dominated 
era after reheating in a model of IV-charged scalars coupled to an abelian gauge field. The spectrum of the magnetic 
field generated during the non-equilibrium evolution was computed using the formulation recently introduced in 
ref. |25l . The dissipative effects of the conductivity are included by separating the contribution from hard modes (with 
momenta of order T ) to the polarization tensor of the gauge fields. These modes are always in local thermodynamic 
equilibrium. The non-perturbative, non-equilibrium dynamics of the scalar field after the phase transition was studied 
in the large N limit. The dynamics after the phase transition features two distinct stages: an early and intermediate 
time, spinodal stage, which is dominated by the growth of long-wavelength fluctuations, followed by a scaling regime 
during which the scalar field becomes correlated over horizon-sized domains. During both regimes, strong non¬ 
equilibrium fluctuations lead to large current-current correlation functions which entail the generation of magnetic 
fields. The scaling regime is the most effective for primordial magnetogenesis since this stage lasts the longest. During 
this stage magnetic fields with superhorizon wavelengths are generated via loop effects, the scalar field momenta in 
the polarization loop corresponds to wavelengths of the order of or shorter than the horizon. Thus causal scalar 
field fluctuations lead to the generation of magnetic fields on superhorizon scales. The generation of magnetic field is 
hindered by the large conductivity of the plasma and equipartition between electric and magnetic fields does not hold. 
The spectrum of the primordial magnetic field is insensitive to the magnetic diffusion length which is sub-horizon 
during the radiation era. 

Our final result for the spectrum generated during the scaling regime is given by eq. (TQ2I) . The ratio of the energy 
density of the magnetic fields on scales larger than L to the energy density in the cosmic background radiation 
r(L,r/) = pb(L,t])/ p cm b(L,r]) is given by equation (14.341) . For values of N, and the gauge coupling consistent with 
particle physics models we find that 


— 34 ( L \ 5 { T EW \ 

^lMpcy \T{rf)J 


r(T(r]),L) ~ 10 


7 


(5.1) 





Therefore, 


24 


r(T( V ),L) 


10 34 at the EW transition 
10 -14 at the QCD transition . 


(5.2) 


Therefore, the large scale primordial magnetic fields generated during the scaling stage after a phase transition may 
be a plausible mechanism to generate primordial magnetic fields which will be further amplified by the collapse of 
protogalaxies and by astrophysical dynamos. 

Probably a phase transition at a temperature much larger than the electroweak leading to a scaling regime lasting 
until the QCD phase transition is ruled out by the temperature inhomogeneities at the last scattering surface. Further¬ 
more the generation of electromagnetic fields on sub-horizon scales requires a full kinetic equation that incorporates 
the microscopic causal processes that lead to thermalization, the study of these is beyond the scope of this article. 

Magnetogenesis after the QCD phase transition? The model that we studied here is assumed to describe 
the robust features from the non-equilibrium dynamics of a charged sector coupled to a (hyper) charge gauge field. 
GUT’s or SUSY theories may provide the corresponding framework. 

However we now argue that precisely the model studied here can actually describe the non-equilibrium dynamics 
after the QCD phase transition(s). After hadronization and chiral symmetry breaking most of the hadrons produced 
will be pions, at least this is the experimental situation in ultrarelativistic heavy ion collisions. Neglecting the charge 
form factor (which is justified for momenta much smaller than the p meson mass rn p ~ 770 Mev) the charged 
pions couple to the electromagnetic held with minimal coupling. The chiral transition is conjectured to be in the same 
universality class as the 0(4) linear sigma model lsfill . Thus the model presented in this article, is the low energy effective 
field theory for the triplet of pions, two charged and one neutral. Thus we conjecture that the study in this article does 
describe the generation of magnetic fields by long-wavelength pions. Therefore the analysis of this article can apply 
to magnetogenesis during the chiral phase transition in QCD. While the charged pions couple to electromagnetism 
via the minimal coupling in the long-wavelength limit, the neutral pion couples to the electromagnetic field through 
the chiral anomaly 7r° —> 2y and such process will also produce magnetic fields. We will study the possibility of large 
scale primordial magnetogenesis during the chiral phase transition in QCD in a forthcoming article. 
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